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Abstract. Supernova 1998bu in the galaxy M96 was observed by COMPTEL for a total of 88 days starting 17 
days after the explosion. We searched for a signal in the 847 keV and 1238 keV lines of radioactive 56 Co from 
this type la supernova. Using several different analysis methods, we did not detect SN1998bu. Our measurements 
should have been sensitive enough to detect 60 Co gamma-rays as predicted from supernova models. Our 2a flux 
limit is 2.3 • 10 -5 photons cm -2 s _1 ; this would correspond to 0.35 Mq of ejected 56 Ni, if SN1998bu were at a 
distance of 11.3 Mpc and transparent to MeV gamma rays for the period of our measurements. We discuss our 
measurements in the context of common supernova models, and conclude disfavoring a supernova event with large 
mixing and major parts of the freshly-generated radioactivity in outer layers. 
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1. Introduction 

Despite their widespread use as 'standard candles', the 
physical nature of supernovae of type la is still not un- 
derstood in terms of physical processes; therefore correc- 



tions of evolutionary aspects remain empirical (Branch 
1998; Niemeyer & Truran , 200C). Type la supernovae 



are believed to be caused by thermonuclear explosions of 



models of the bolometric light curve of type la supernovae 
imply that typically 0.3-0.5 M Q of ra dioactive 56 Ni energy 
are needed ( Hoflich fc Khoklov, 1996| ). Observations of the 
supernova light curve's peak magnitude or of NIR lines 
of Fe[II] and Co[H] have mostly been used to derive 56 Ni 
masses; the rather wide range of inferred 56 Ni (0.1 - 1.14 
M0)(Contardo et al., 2001) is difficult to understand if a 



CO white dwarfs flLiyio, 2000J [Hoflich fc Khoklov, 1996| ; 



Nomoto et al., 1997). Radioactive energy of ~0.5 M 



of 56 Ni synthesized in such explosions is considered the 
driver of all types of observed light from these objects. 
The nearby SN1998bu supernova was a unique opportu- 
nity to directly measure gamma-rays from the 56 Ni decay 
chain with the Compton Observatory. 

The dynamics of a white dwarf explosion are difficult 
to model due to the ran ge of scales involved i n flame ig- 
nition and propagation (Iwamoto et al., 1999). Theories 
which are most successful in describing observations of 
type la supernovae include empirical components for key 
aspects. Observationally type la supernovae are a fairly 
homogeneous phenomenon ( Branch, 1998 ), suggesting a 
narror range of synthesized 5ti Ni masses. Constraints from 
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single well-tuned process is held responsible for the super- 
novae of type la and in particular their " standard candle" 
characteristics. Among different types of models, a rather 
wide range of 56 Ni masses of 0.1 up to 0.8 M Q is discussed 
(see 'Discussion' section below). 

Critical parameters of the explosion models are the ig- 
nition density of the white dwarf at its core, and the tran- 
sition from the early deflagration stage (sub-sonic flame 
propagation) into a detonation (super-sonic flame prop- 
agation)^. g. Lcibundgut, 2000; Livio, 2000). The former 
is estimated from evolutionary models of the (uncertain) 
progenitor and the effective binary accretion rate, involv- 
ing uncertain issues such as steady or flash-like nuclear 
burning of the accreted H and He material, or modula- 
tion of the accretion flow from the companion through 
the wind of the white dwarf (at the higher accretion rates 
required by the lower ignition densities preferred from oth- 
erwise excessive production of neutron-rich Fe group iso- 
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topes(Nomoto ct al., 1997)). The acceleration of the flame 
speed can only be treated as a purely empirical param- 
eter of models at present, but critically determines the 
final 56 Ni mass and the Fe group to lighter element ratio 
(Iwamoto et al., 1999). Three-dimensional model treat- 
ments of this flame "micro-physics" is promising, but still 



a challenging problem (Reinecke et al., 1999). 

The fact that type la supernovae are rather homoge- 



neous (Branch, 1998) suggests a clear evolutionary path 
towards a well-definied presupernova star and a robust ig- 
nition condition. This led to the model of binary accretion 
of H and He rich matter onto a CO white dwarf at a well- 
tuned accretion rate such that H and especially He nu- 
clear burning proceeds non-catastrophically and the white 
dwarf C mantle grows in mass until the Chandrasekhar 
mass limit is reached; the thermonuclear runaway explo- 
sion ensues from fast nuclear burning of Carbon ignited 
at the core due to heating from gravitational pressure and 
H and He shell-burning heat conduction ( Nomoto, 1982|). 



and its evolution reflect the supernova envelope structure, 
with much less information on the core. Spectral informa- 
tion tells us about material mixing and the total kinetic 
energy. However optical photons are created long after the 
initial explosion; most information from the early stage 
of the supernova event is lost. This makes it difficult to 
discriminate between different explosion models or model 
parameters. 

Observations of 7-rays promise more direct informa- 
tion from the core of the supernova and the explosion 
mechanism. The radioactivity produced in the initial 
event decays and produces gamma-ray lines, which can 
be observed directly once the supernova is transparent 
to gamma-rays (after about 30-100 days). Even in the 
early stages, 7-rays from radioactivity will escape from 
outer layers, their intensity depending on ejecta mixing. 
Differences in the predicted 7-ray spectra have been sug- 
gested as the key observation to discriminate between 
models and the extent of mixing (Hoffich et al., 1998; Pintc 



Sub-Chandrasekhar mass white dwarfs could also explode et al., 2001). This is best observed at early times; at times 



as type la supernovae ( Livne, 1990] ): The merging of two 



white dwarfs would disrupt the lighter of the two into a C- 
O envelope accreting onto the more massive white dwarf 
("double-degenerate" model). As the accretion proceeds 
to exceed the Chandrasekhar limit, the white dwarf ig- 
nites C centrally as above (Iben & Tutukov, 1984). There 
are some doubts if the merging process will avoid core 
collapse of the merged object and produce a thermonu- 
clear supernova; e.g. transport of rotational energy is crit- 
ical (Livio, 200C). Alternatively, a single-degenerate sub- 
Chandrasekhar model has been proposed: A He layer built 
up from accretion and steady hydrogen burning as above 
may ignite in a flash and thus send a shock wave into 
the white dwarf core, adding to the gravitational heat 
and thus also ignit ing carbon in the center for a lower- 



later than 30 days, differences in total 56 Ni masses can 
mimic differences between sub- and Chandrasekhar mod- 



mass white dwarf (|Nomoto, 1982| |Livnc, 1995| ). In this 



scenario much of the radioactive material would be pro- 
duced towards the outside, resulting in different evolu- 
tion of radioactivity-derived supernova light. Recent con- 
straints on early spectra and the absence of intermediate- 
mass elements in the outer fast ejecta disfavor this scenario 



somew lat ( Livio, 2000 )_. 

A deflagration model was recently favored for 



els ( |Pinto et al., 200l| ). 

The sensitivity of the 7-ray instruments on-board 
CGRO (OSSE and COMPTEL) limits the observations 
of type la supernovae to events within about 15 Mpc. To 
date, only one ev ent, SN1991T, was marginally detected 
with COMPTEL ( |Morris et al., 1997[ ). SN1998bu provides 
a second opportunity for line searches in type la super- 
novae. Some theoretical models predict 7-ray line fluxes 
well above the sensitivity limits of COMPTEL and OSSE 
for an assumed distance of 11 Mpc. Although COMPTEL 
lacks the spectral resolution to provide unique and deci- 
sive 7-ray line shape diagnostics, an independent proof of 
the radioactive 56 Ni mass origin through detection of the 
corresponding 7-ray line fluxes was attempted, and is im- 
portant given the complexity and unknowns of conversion 
of radioactive energy in a supernova envelope. 

On May 9.9 UT in 1998 supernova S N1998bu w as dis- 
covered in the galaxy M96 (NGC 3368) flVilli, 1998[ ) From 
wide-band spectrograms it was classified as type la ( Ayani 
ct al., 1998|; Meikle et al., 1998 ) . From an earlier observa 



SN1998bu on purely spectroscopic arguments (Vinko 
et al., [200 1| ). And the occurrence near or in a spiral arm 



and the observation of a light-echo ( Cappellaro et al 
200l|)~}nay suggest that the progenitor system could be 
younger than those of the average type la supernovae, 
therefore anomalous and igniting at a particularly low 
density. 

With such diversity of models and the difficulties of 



tion ( Faranda fe Skiff, 199S ) and an estimate of the max- 
imu m blue light at £ B m ax = 10952.7 ± 0.5 TJD (i.e. May 
19), Meikle et al. (1998) inferred the date of the explosion 
to be May 2.0 ± 1.0 UT (i.e. TJD 10935 ± 1). The distance 
to M96 had been known from HST Ccphcid measurements 
as 11.3 ± 0.9 Mpc QHjorth fe Tanvir, 1997| ), though an- 
other value of 9.6 ± 0.6 Mpc had been derived from mea- 



surements of planetary nebulae (Feldmeier et al., 1997) 



SN1998bu appears to be a typical type la event ( Jah et al 



detailei physics modeling, observations of a variety of as- 1999), its reddening can be attributed to dust in the host 



pects of type la supernovae are a key to clarify the true 
nature of these events. 

Observations of a large sample of supernovae in UV, 
optical and infrared bands have been made and discussed 
widely. But this radiation originates from driving pro- 
cesses deep inside the object, the bolometric light curve 



galaxy. From its dereddened bri ghtness it was concl uded 
(using the methods d escribed in Nomoto ct al. (1997 ) and 
Iwamoto ct al. (1999 )) that the total ejected Ni mass was 
rather typical. A value of 0.77 M Q has been derived from 



analysis of the bolom etric light curve ( Lcibundgut, 200 
IContardo et al., 200 1| ). 
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Fig. 1. The total energy spectrum (left) and the background subtracted spectrum (right) for a sample position on the 
grid around SN1998bu. 
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Fig. 2. The result of the spectral analy- 
sis. The derived 2a upper limits are 4.1- 
1 000 10" 5 photons cm" 2 s" 1 for the 847 keV 
line and 2.3 TO -5 photons cm -2 s" 1 for 
the 1238 keV line. 



2. Observations 

The COMPTEL observations of SN1998bu began on TJD 
10952, 17 days after the explosion. Due to this late start 
and to the low sensitivity of COMPTEL for low energies, 
we missed the opportunity to measure the decay lines 
of 56 Ni (r = 8.8 d) at 750 keV and 812 keV. 56 Co, the 
daughter nucleus of 56 Ni, has a mean life of r = 112 d 
and decays to stable 56 Fe. In this decay two 7-ray lines 
are emitted at 847 keV and 1238 keV. The COMPTEL 
observations spanned a total duration of 88 days, ending 
at TJD 11071, 136 days after the explosion. 

The COMPTEL telescope detects 7-rays through a 
Compton scatter interaction in its upper plane of liquid- 
scintillation detectors ("Dl"), followed by detection of an- 
other interaction in the lower plane of Nal(Tl) scintillation 
detectors ("D2"), ideally totally absorbing the scattered 
7-ray. Imaging information is provided by the Compton 



scatter kinematics through energy deposit and interac- 
tion location measurements (For a detailed description of 
COMPTEL see |Schonfelder et ah, 1993|). 



3. Analysis and Results 

To obtain a higher sensitivity at lower energies the ob- 
servations were carried out in a "low-threshold" mode, 
where the energy thresholds of the D2 modules were re- 
duced below the normal 600 - 700 keV range; the normal 
mode would have degraded the instrument response for 
the 847 keV photons. Two different analysis approaches 
were used: "spectral" and "imaging" analysis. 

3.1. Spectral analysis 

In this approach the full spectral response, i.e. the photo 
peak and the Compton tail, is exploited. Imaging param- 
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eters of the measurement are used for coarse field-of-view 
selection of the events. Only events with D2 energy de- 
posits above 600 keV were used to avoid the 511 keV line 
background which originates from positrons. These result 
mainly from 22 Na within the instrument, one of the major 
sources of instrumental background, which has its origin 
in activation of COMPTEL passing through the radiation 
belts (Weidcnspointner et al., 2001). 



Spectra were derived for a ~ 12° FWHM beam on a 
grid around the SN position with a grid spacing of 3°. 
Fig. [l] shows an example for the direction covering the 
SN position. The beam size was chosen through the se- 
lections of events from within a 2° wide range around 
the pivot direction, performed near the cone-shaped re- 
sponse in the COMPTEL data space. Background spectra 
were constructed from the same direction, selecting a ring- 
like 3° — 7° wide region around the same pivot direction. 
Background-subtracted spectra were then analyzed for the 
relevant Co decay lines (see Fig. |l|). Gaussian-shaped lines 
with instrumental resolution are fitted to the data to de- 
termine line intensities. This analysis and in particular the 
background determination may include unknown system- 
atic effects that could mimic lines. To account for these un- 
certainties (in addition to the statistical uncertainties) , we 
determined the variance of the fitted line intensities across 
our measured sample empirically: We histogrammed the 
fitted line intensities from all spatial grid points and used 
the width of this distribution to estimate uncertainties and 
significances, for both 7-ray lines of interest (see Fig. |f). 
The exposure and effective detector area were used to con- 
vert line intensities and uncertainties into source photon 
flux units. 

No significant difference between the spectra of the 
supernova and off-supernova positions was found. We ob- 
tained a 2<T upper flux limit of 4.1-10 -5 photons cm~ 2 s _1 
for the 847 keV line and 2.3 • 10~ 5 photons cm" 
the 1238 keV line. 



3.2. Imaging analysis 

In this analysis the imaging capability of COMPTEL was 
fully exploited, thus increasing the sensitivity of the in- 
strument mainly in the 847 keV line due to a better 
background suppression. In the "low-threshold" observing 
mode most of the D2 detector thresholds are well below 
the 847 keV line of interest. However, there is a consider- 
able spread of thresholds, the lowest value being 450 keV. 
In standard COMPTEL analysis, event selection with D2 
energy deposits above 650 keV is applied, to simplify the 
modular composition of the instrument response function. 
This is now inappropriate, and response function composi- 
tion and usage is more complex, as a price for the improved 
847 keV sensitivity. For the upper detector plane Dl, all 
modules have hardware thresholds in a narrow range at or 
below 50 keV, thus still allowing the use of a homogeneous 
analysis threshold at 50 keV for all Dl detectors. 



Instrument imaging response functions were derived 
for our analysis from simulations, where each D2 module 
was treated with its effective hardware threshold. The full- 
instrument data were analyzed in two different ways: We 
first used a summed response function for both lines, but 
for the 847 keV line we also analyzed separately for "mini- 
telescopes" and composed the individual results. 

We demonstrate the difference between standard anal- 
ysis and our more complex "low-energy-threshold-mode" 
response function on 3 days of Crab low-threshold data, 
fortuitously collected during an observation of Geminga. 
Using the low-mode response we clearly detect Crab with 
4.5(7 at 847 keV, compared to only 3. Oct with the stan- 
dard response (see Fig. ||). For the 1238 keV line the val- 
ues are 1.6ct and 0.63ct, respectively. This shows that most 
D2 modules are now sensitive below the 650 keV analy- 
sis threshold used for the standard analysis, resulting in a 
larger sensitivity mainly for the 847 keV line. 

Unfortunately, due to the strong background line at 
511 keV, we did not obtain the theoretically- achievable 
sensitivity at the 847 keV line. To reduce the impact of 
this background via its scattering pattern, we applied a cut 
of 40° in the angle between the incident and the scatter ed 
7-ray (the so-called (p angle (see Schonfelder et al., 1993| )). 
This sacrifices a fraction of source events. 

Using these simulated summed responses and the (p 
cut, maps for the two 7-ray energies (shown on the left side 
in Figure ^J) are produced through a maximum-likelihood 
method. There is no detectable signal from the supernova. 

For determination of an upper limit on the flux in the 
7-ray lines we estimated statistical and systematic uncer- 
tainties. In a histogram of flux values for all pixels in the 
maps of Fig. ||, the flux variances provide a global measure 
of the uncertainty. Using the Bayesian method described 
in Georgii et al. (1997), which accounts for the systematic 
and statistical uncertainties, 2ct upper limits of 5.5 • 10 -5 
photons cm~ 2 s _1 for the 847 keV line and of 3.2 • 10 -5 
photons cm~ 2 s _1 for the 1238 keV line were determined. 
The less constraining limit for the 847 keV line is a result 
of the tp cut for background suppression. 

Therefore another analysis method was applied to the 
847 keV line, only. Here separate maximum-likelihood 
imaging of the SN1998bu region was made for each mini- 
telescope (the combination of one single D2 module with 
all Dl modules) by applying the specific response for each 
D2 module. With energy cuts at 650 keV, but no (p cut, 
511 keV background suppression retains better 847 keV 
sensitivity than above. Since the likelihood values are ad- 
ditive (probabilities) we combined these results to derive a 
2ct upper limit of 3.1 • 10~ 5 photons cm~ 2 s _1 (see Fig. ||). 



3.3. Summarizing the results 

From the imaging analysis we obtain a 2ct upper limit 
of 3.1 • 10~ 5 photons cm~ 2 s _1 for the 847 keV line and 
3.2 -10~ 5 photons cm" 2 s" 1 for the 1238 keV line. Spectral 
analysis yields 4.1 • 10~ 5 photons cm~ 2 s _1 for the 847 keV 
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Fig. 3. Improvement of the sensitivity with the low-threshold mode response function (left) for 3 day of Crab data 
compared to the sensitivity with the standard response function (right). Shown is a significance map based on a 
maximum-likelihood method, a contour lines are plotted. The position of Crab is marked with a star. 
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Fig. 4. Flux maps of the SN1998bu region in 847 keV (top left) and 1238 keV (bottom left). No significant excess is seen 
around the position of the SN, which is marked with a star. The contour levels are in units of 1CP 6 photons cm~ 2 s — 1 . 
Right to each map its flux distribution and a Gaussian fit to it is shown. The distortion of the distributions around 
small flux values is due to the method applied. 



line and 2.3 • 10 -5 photons cm -2 s _1 for the 1238 keV line, 
respectively. The reason for the different values for both 
lines in the spectral analysis lies in the steep decrease of 
the background with the energy. In contrast, in imaging 
analysis the background on the sky around the source is 
not much different for both lines. 



The analysis presented here was based on the assump- 
tion that the emitted lines are narrow. This is justified 
because of the energy resolution of COMPTEL (about 88 
keV for the 847 keV and about 110 keV for the 1238 keV 
line). However, due to the velocity of the expelled matter, 
the lines are broadened and even shifted (only shortly af- 
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Fig. 5. The flux limit for the 847 keV 
line, derived from combined imaging re- 
sults per each minitelescopes. We fix 
a source flux at the position of the 
SN, and determine the maximum log- 
likelihood ratio value for each detector 
subset. We can add these, as derived 
from independent data. We normalize 
to the log-likelihood ratio of our best- 
fit flux value of photons cm -2 s _1 . 
Varying adopted fluxes, we derive the 
parabolic behavior near the minimum 
(see stars). For our 11 independently- 
varied parameters (one per each detec- 
tor subset) we obtain a 2a SN flux 
limit by increasing the likelihood by 
22 units above the minimum (i.e., 3.1 • 
10 _5 photons cm -2 s _1 ). 
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Fig. 6. The model fluxes for the 847 keV (left) and for the 1238 keV (right) line for different models versus time after 
the explosion for a distance of 11.3 Mpc. The upper limits from the spectral and imaging analysis are also shown. 
The solid line represents the more sensitive value from both methods for each line. Note that for the 847 keV line the 
imaging analysis and for the 1238 keV line the spectral analysis is more sensitive. 



ter the explosion) by the Doppler effect (see Hofiich et al 
(1998|)j |Iscrn (1997 )). An analysis of the expected degra- 
dation of the upper limits yields for a persistent Doppler 
broade ning of the 1238 keV line of 10 000 km/sec flHoflich 
et al., f 998 ) , an increase of the 2a upper limit flux by 12 
% compared to the narrow-line assumption. This changes 
the flux upper limits only marginally and we therefore ne- 
glected it. 

4. Discussion 

To compare 7-ray upper limits with theoretical Ni mass 
predictions, the distance to the SN plays an essential role. 
The host galaxy M96 ha d a HST-Cepheid-d etermined dis- 
tance of 11.6 ± 0.9 M pc (Tanvir et al., 1995 ), later revised 
to 11.3 ± 0.9 Mpc by [Hjorth fc Tanvir (1991 ). This makes 
SN1998bu one of seven SN observed in galaxies with a 



distance set by Cepheid measurements. We adopt this dis- 
tance of 11.3 Mpc for our analysis. 

A second key factor is the transparency of the su- 
pernova to 7-rays. This is a key issue determining how 
much radioactive energy is converted into kinetic energy 
and supernova light. The maximum of the optical light 
curve (about 10 days after the explosion) is set by a max- 
imum of the product of (declining) energy deposition and 
(rising) 7-ray energy escape (Pinto and Eastman, 2001). 
Furthermore, the Compton scattering optical depth to 1 
McV 7-rays is below unity beyond about 50 days after 



explosion (Pinto et al., 2001), but absorption corrections 
to observed line 7-rays are probably significant for typical 
mode ls up to about 100 d ays after explosion (see e.g. Fig. 
11 in Hofiich et al., 1998, calculated for energies down to 



10 keV, however). This illustrates clearly the importance 
of 7-ray measurements with high spectral resolution and 
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at those early times, in order to directly address the ex- 
plosion mechanism: the line shapes and the ratio of the 
different line intensities from the 56 Ni decay chain can 
reveal the ratio between deposited and directly radiated 
radioactive energy (e.g. [Hoflich fc Khoklov (1996; )). 

For illustrative purposes and simplification, we may 
simply assume as an extreme case that the supernova 
was transparent for our observation of the gamma-rays 
from 56 Co decay over days 17-136 with emphasis on the 
late part; in this case we directly convert our flux lim- 
its into 56 Co (and therefore original 56 Ni) masses. Our 
lowest upper limit for the 1238 keV line of 2.3 • 10 -5 
photons cm~ 2 s _1 then constrains the visible 56 Ni mass 
to below 0.35 M Q . If we then want to reconcile this 
with the 0.77 M o f total 56 Ni determined bolometrically 
( |Leibundgut, 200(]| ), more than half of the 7-ray energy 
would be deposited in the supernova over this time win- 
dow. We therefore do have to look in detail at the en- 
ergy deposition efficiency around peak optical luminosity 
and/or effectiveness of 7-ray escape soon thereafter. 

For several model classes (detonation, delayed deto- 
nation, and sub-Chandrasekhar), 7-ray light curves have 
been calculated through detailed Monte-Carlo photon 



transport in the expanding supernova (Hoflich et al 



199S; kumagai, 1998; tsern, 1997; Pinto et al, 2001) 



Considerable variety in the gamma-ray flux by factors up 
to 5 arises from the different explosion models, envelope 
structures, and photon transport treatments employed in 
such calculations. In Fig . p|, expecte d 7-ray light curves 
for a few typical models (Isern, 1995; Kumagai, 1998[)_ are 
shown, re-scaled for a distance of 11.3 Mpc. In Table [l] we 
list the 56 Ni mass for each of these models, together with 
time-averaged fluxes over the observation time for each 7- 
ray line. We see that the predicted 56 Co 7-ray flux does 
not follow the straightforward scaling to the amount of 
56 Ni, the explosion mechanism and the envelope photon 
transport determine the time-dependent 7-ray fluxes. For 
the same type of explosion model, predicted 56 Ni masses 
vary within a factor of two: For the delayed-detonation 
class of models, values between 0.55 M s and 0.96 M 



have b een published ( Iwamoto et al., 199E ; Wooslcy and 



Weave: -, 1991 ; tsern, 1997 ) , as a result of differences in the 
point at which the initially-slow nuclear burning (defla- 



gratio n) is assumed to turn into a detonation. This typi- 



cal intrinsic variability within an explosion type of a fac- 
tor of two, which directly translates into the 7-ray flux 
scaling, indicates the systematics which typically remains, 
within an explosion type. In fact, any of the SNe la sce- 
narios (sub-Chandrasekhar, deflagration, delayed detona- 
tions, and pulsating delayed detonation models) has been 
shown to be capabable to produce a wide variety of 56 Ni 
masses ranging from ~ 0.1 to 1 M 



I (e.g. 


Nomoto et al. 


H6fiich et al. (199?)). 



On the other hand, Fig. || illustrates that, for approx- 
imately the same amount of total 56 Ni, pure deflagration 
or detonation models are about a factor of two dimmer in 
7-rays, while the sub-Chandrasekhar model with substan- 



tial 56 Ni sitting further outside reaches a 7-ray flux about 
twice as large as the typical delayed-detonation model. 

Table 1. Fluxes in the 847 keV and 1238 keV line for 
different models, averaged over our observation time. 



Model 


bb Ni mass 


Is47 keV-10 


Il238 kcv40 b 




[M s ] 


[ph cm _2 s _1 ] 


[ph cm~ 2 s _1 ] 


W7f] 


0.58 


4.2 


3.0 


W7DT a 


0.77 


5.8 


4.1 


HeCD a 


0.72 


8.2 


5.5 


WDD2 a 


0.58 


3.9 


2.8 


Deflag.| 


0.50 


1.5 


1.0 


Del. Det. b 


0.80 


4.3 


2.7 


Det. 6 


0.70 


4.0 


2.7 


SubCH. 6 


0.60 


2.7 


1.5 



Kumagai (199^ ) 



tsern (1998) 



In Table [l] we also list the time-averaged 56 Co 7-ray 
fluxes of each model, together with our 2er upper limits. 
Our SN1998bu flux limits are well below the "HeCD" and 
the "W7DT" model predictions for both lines. The "W7", 
" WDD2" , delayed detonation and detonation model fluxes 
are marginally consistent with our flux limits, while the 
fluxes predicted from the sub-Chandrasekhar model and 
the deflagration model are consistent with our limits for 
both line energies at the adopted distance. 

This comparison illustrates that with a 56 Ni mass in 
the " typical" range der ived for SN1998bu, around 0.7-0.8 
Mq ( Leibundgut, 2000 ), we should have seen 56 Co 7-rays 
at least due to those models which turn more rapidly from 
deflagration into detonation (W7DT) or partially-produce 
radioactivity in their outer ejecta (HECD). Yet, within 
Chandrasekhar-type models of the presently-favored type 
of a delayed transition from deflagration into detonation, 
a total 56 Ni mass as high as about 1 M may still be 
consistent with our measurement. 



At a distance of 9.6 ± 0.6 Mpc ( [Fcldmcicr et al., 1997p 
based on planetary nebulae (PN) all models would be in- 
consistent with our 1238 keV and 847 keV flux limits, in- 



dicat ing either this distance is incorrect (see ( Maoz et al 



1999 ) for a discussion on a possible correction in the dis- 
tance ladder scale) or that model treatments generally 
overestimate the 56 Ni masses. 

It will require time-resolved measurements of the 7-ray 
flux (hence a brighter / more nearby supernova or a more 
sensitive instrument), or exploitation of spectral-shape de- 
tails as promised by the spect rometer aboard INTEGRAL 
(see discussion in [scrn, 1997 ), to decide among explosion 
models from gamma-ray line measurements alone. 
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